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Abstract 

Quasi-realistic string models in the free fermionic formulation typically 
contain an anomalous U{\), which gives rise to a Fayet-Iliopoulos Z)-term 
that breaks supersymmetry at the one-loop level in string perturbation theory. 
Supersymmetry is traditionally restored by imposing F~ and D-flatness on the 
vacuum. By employing the standard analysis of flat directions we present a 
quasi-realistic three generation string model in which stringent F- and D- 
flat solution do not appear to exist to all orders in the superpotential. We 
speculate that this result is indicative of the non-existence of supersymmetric 
flat F- and D-solutions in this model. We provide some arguments in support 
of this scenario and discuss its potential implications. Bose-Fermi degeneracy 
of the string spectrum implies that the one-loop partition function and hence 
the one-loop cosmological constant vanishes in the model. If our assertion 
is correct, this model may represent the first known example with vanishing 
cosmological constant and perturbatively broken supersymmetry. We discuss 
the distinctive properties of the internal free fermion boundary conditions that 
may correspond to a large set of models that share these properties. The 
geometrical moduli in this class of models are fixed due to asymmetric boundary 
conditions, whereas absence of supersymmetric flat directions would imply that 
the supersymmetric moduli are fixed as well and the dilaton may be fixed by 
hidden sector nonperturbative effects. 
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1 Introduction 



The quasi-realistic heterotic-string models in the free fermionic formulation, 
which are related to Z2 x Z2 orbifold compactifications, are among the most re- 
alistic string models constructed to date. These models produce a rich variety of 
three generation models with the canonical 5*0(10) embedding of the Standard Model 
spectrum, and include: the flipped SU{5) string models [1] (FSU5); the standard-like 
string models [2, 3, 4, 5, 6, 7]; the Pati-Salam string models [8] (PS); the Left-Right 
symmetric string models [9] (LRS). Many of the issues pertaining to the phenomenol- 
ogy of the Standard Model and Grand Unification were investigated in the context of 
these models [10]. Additionally, the free fermionic models produced the first known 
string models in which the matter content in the observable Standard Model charged 
sector of the effective low energy quantum field theory consists solely of the Minimal 
Supersymmetric Standard Model [5, 6]. 

A common feature of many of the quasi-realistic free fermionic heterotic-string 
models is the existence of an "anomalous" U{1) symmetry [11]. The anomalous U{1)a 
is broken by the Green-Schwarz-Dine-Seiberg-Witten mechanism [12] in which a 
potentially large Fayet-Ihopoulos D-term ^ is generated by the VEV of the dilaton 
field. Such a D-term would, in general, break supersymmetry, unless there is a 
direction = J2oii4>i in the scalar potential for which I]Q^|Q;ip is of opposite sign 
to ^ and that is D-fiat with respect to all the non-anomalous gauge symmetries, 
as well as F-flat. If such a direction exists, it will acquire a VEV, cancelling the 
Fayet-Iliopoulos ^-term, restoring supersymmetry and stabilising the vacuum. The 
set of D- and F-flat constraints is given by 
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{Da) = {D^) = 0; {Fi = -^) = ; (1.1) 
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where Xk are the fields which acquire VEVs of order while the i^'-terms contain 
fields Tji like squarks, sleptons and Higgs bosons whose VEVs vanish at this scale. 
Q\ and denote the anomalous and non-anomalous charges, and Mpi fa 2 x 10^^ 
GeV denotes the reduced Planck mass. The solution [i.e. the choice of fields with 
non-vanishing VEVs) to the set of equations (1.1)-(1.3), though nontrivial, is not 
unique. Therefore in a typical model there exist a moduli space of solutions to the 
F and D fiatness constraints, which are supersymmetric and degenerate in energy 
[13]. Much of the study of the superstring models phenomenology (as well as non- 
string supersymmetric models [14]) involves the analysis and classification of these 



2 



flat directions. The methods for this analysis in string models have been systematised 
in [15, 16, 5, 17]. 

In general it has been assumed in the past that in a given string model there should 
exist a supersymmetric solution to the F and D flatness constraints. The simpler type 
of solutions utilise only fields that are singlets of all the non-Abehan groups in a given 
model (type I solutions). More involved solutions (type II solutions), that utilise also 
non-Abelian fields, have also been considered [17], as well as inclusion of non-Abelian 
fields in systematic methods of analysis [17]. The general expectation that a given 
model admits a supersymmetric solution arises from analysis of supersymmetric point 
quantum field theories. In these cases it is known that if supersymmetry is preserved 
at the classical level, i.e. tree-level in perturbation theory, then there exist index 
theorems that forbid supersymmetry breaking at the perturbative quantum level 
[18]. Therefore in point quantum field theories supersymmetry breaking may only be 
induced by non-perturbative effects [19]. 

Recently we constructed string models [7] in which the issues of supersymmetric 
flat directions merits further investigation. The aim of ref. [7] was to build models 
with a reduced untwisted Higgs spectrum. This was achieved in ref. [7] by impos- 
ing asymmetric boundary conditions in a boundary condition basis vector that does 
not break the 5*0(10) symmetry. An unforeseen consequence of the Higgs reduction 
mechanism of ref. [7] was the simultaneous projection of untwisted SO (10) singlet 
fields. Subsequently the moduli space of supersymmetric flat solutions is vastly re- 
duced. In fact, in ref. [7] it was concluded that the model under investigation there 
does not contain supersymmetric flat directions that do not break some of the Stan- 
dard Model symmetries. Indeed, for that reason, a phenomenologically viable model 
with a reduced untwisted Higgs spectrum was not presented in ref. [7]. 

The question therefore remains whether a phenomenologically viable model with 
a reduced untwisted Higgs spectrum exists. In this paper we explore this question 
further. The untwisted Higgs reduction mechanism that we use here differs from the 
one of ref. [7] . Here we present a model that utilises boundary conditions that are both 
symmetric and asymmetric in the basis vectors that break 5*0(10) to 50(6) x 50(4), 
with respect to two of the twisted sectors of the Z2 x Z2 orbifold. The consequence is 
that two of the untwisted Higgs multiplets, associated with two of the twisted sectors, 
are projected entirely from the massless spectrum. As a result, and similar to the 
model of ref. [7], the string model contains a single pair of untwisted electroweak 
Higgs doublets. 

In the process of seeking such a model with a phenomenologically viable super- 
symmetric flat direction, we arrive in this paper to the unexpected conclusion that 
the model may not contain supersymmetric flat directions at all. In the least, this 
model appears to have no D-flat directions that can be proven to be F-flat to all 
order, other than through order-by-order analysis. That is, there does not appear 
to be any D-Hai directions with stringent F- flatness (as defined in [5, 6, 20]). In 
the analysis of the fiat directions we include all the fields in the string model, i.e. 
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Standard Model singlet states as well as Standard Model charged states. The model 
therefore does not contain a D-flat directions that is also stringently F-flat to all 
order of non-renormalizable terms. 

The model may of course still admit non-stringent flat directions that rely on 
cancellations between superpotential terms. However, past experience suggests that 
non-stringent flat directions can only hold order by order, and are not maintained to 
all orders [21, 9]. This is the key difference between the string theory case, in which 
heavy string modes generate an infinite tower of terms, versus the field theory case in 
which heavy modes are not integrated out. We therefore speculate that in this case 
supersymmetry is not exact, but is in general broken at some order. 

If this finding remains true after the entire parameter space of possible all-order 
non-stringent fiat directions has been examined (an undertaking of several years), 
we must ask what arc the the implications. If a model without all-order F-fiatness 
were to be found, then supersymmetry would remain broken by the Fayet-Iliopoulos 
term at a finite order, which is generated at the one-loop level in string perturbation 
theory, rather than be cancelled by a D-Hat direction with anomalous charge. If so, 
then this would imply, although supersymmetry is unbroken at the classical level and 
the string spectrum is Bose-Fermi degenerate, that supersymmetry may be broken 
at the perturbative quantum level. Nevertheless, since the spectrum is Bose-Fermi 
degenerate, the one-loop cosmological constant still vanishes. 

The string model that we present contains three chiral generations, charged under 
the Standard Model gauge group and with the canonical 5*0(10) embedding of the 
weak-hypcrchargc; one pair of untwisted electroweak Higgs doublets; a cubic level 
top-quark Yukawa coupling. The string model therefore shares some of the phe- 
nomenological characteristics of the quasi-realistic free fermionic string models. It 
may therefore represent an example of a quasi-realistic string model, with vanishing 
one-loop cosmological constant and perturbatively broken supersymmetry. 

Our paper is organised as follows: in section 2 we review some aspects of the free 
fermionic formalism and in section 3 we elaborate on stringent fiat directions. Then in 
section 4 we present the string model. In section 5 we discuss the methodology of flat 
direction analysis and the evidence for concluding that supersymmetric (stringent) 
flat directions may not exist to all order in the string model of section 4. Section 6 
concludes the paper. 

2 Free Fermionic Models 

In this section we briefly review the construction and structure of the free 
fermionic standard like models. The notation and further details of the construction 
of these models are given elsewhere [2, 3, 22, 5, 9, 7]. In the free fermionic formulation 
of the heterotic string in four dimensions [23] all the world-sheet degrees of freedom, 
required to cancel the conformal anomaly, are represented in terms of free fermions 
propagating on the string world-sheet. In the light-cone gauge the world-sheet 
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field content consists of two transverse left- and right-moving space-time coordinate 
bosons, 2 and X121 and their left-moving fermionic superpartners, and addi- 
tional 62 purely internal Majorana-Weyl fermions, of which 18 are left-moving and 
44 are right-moving. The models are constructed by specifying the phases picked by 
the world-sheet fermions when transported along the torus non-contractible loops 

(2.1) 

Each model corresponds to a particular choice of fermion phases consistent with 
modular invariance and is generated by a set of basis vectors describing the transfor- 
mation properties of the 64 world-sheet fermions. The physical spectrum is obtained 
by applying the generalised GSO projections. The low energy effective field theory 
is obtained by S-matrix elements between external states [24] . 

The boundary condition basis defining a typical realistic free fermionic heterotic 
string model is constructed in two stages. The first stage consists of the NAHE 
set, which is a set of five boundary condition basis vectors, {1, 5, 61, 62, ^3} [25, 22]. 
The gauge group, after imposing the GSO projections induced by the NAHE set, is 
5'O(10) X SO{e>f X E^, with = 1 supersymmetry. The NAHE set divides the inter- 
nal world-sheet fermions in the following way: generate the hidden i^g gauge 
group, generate the 5*0(10) gauge group, while {yT"''"'^ ^f]^}-, {y^, y^, cj^, cj®, 77^} 
and {uj^'"''^,ff} generate the three horizontal 5*0(6) symmetries. The left-moving 
{y.,uj} states are divided to {y^''"'^}, {y^,y'^,LO^ ,lo^}, {uj^-'"-'^}, while x^^, x"^-, X"^ 
generate the left-moving N — 2 world-sheet supersymmetry. 

The second stage of the basis construction consists of adding to the NAHE set 
three additional boundary condition basis vectors. These additional basis vectors 
reduce the number of generations to three chiral generations, one from each of the 
sectors bi, 62 and 63, and simultaneously break the four dimensional gauge group. 
The assignment of boundary conditions to {{jj^''"'^} breaks 5'O(10) to one of its sub- 
groups. Similarly, the hidden Eg symmetry is broken to one of its subgroups. The 
flavour 5*0(6)'^ symmetries in the NAHE-bascd models arc always broken to flavour 
U (1) symmetries, as the breaking of these symmetries is correlated with the number 
of chiral generations. Three such U{l)j symmetries are always obtained in the NAHE 
based free fermionic models from the subgroup of the observable Es, which is orthog- 
onal to 5'O(10). These arc produced by the world-sheet currents fi^f]^ (j = 1,2,3), 
which are part of the Cartan sub-algebra of the observable Eg. Additional unbro- 
ken U{1) symmetries, denoted typically by U{l)j (j = 4,5, ...), arise by pairing two 
real fermions from the sets {y^' "'^}, {y^'^, c)^'®} and {Cu^' " ''^}. The final observable 
gauge group depends on the number of such pairings. Alternatively, a left-moving 
real fermion from the sets {y^''"'^}, {y^'"^ , u^'^} and {u^''"'^} may be paired with its 
respective right-moving real fermion to form an Ising model operator, in which case 
the rank of the right-moving gauge group is reduced by one. The reduction of un- 
twisted electroweak Higgs doublets crucially depends on the pairings of the left- and 
right-moving fermions from the set {y, o)}^"""^. 
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Subsequent to constructing the basis vectors and extracting the massless spec- 
trum, the analysis of the free fermionic models proceeds by calculating the super- 
potential. The cubic and higher-order terms in the superpotential are obtained by 
evaluating the correlators 

Am^{vMv^---V'n), (2.2) 

where iV^) are the fermionic (scalar) components of the vertex operators, using 
the rules given in [24]. Typically, one of the C/(l) factors in the free-fermion models 
is anomalous and generates a Fayet-Iliopoulos term which breaks supersymmetry at 
the Planck scale [12]. A supersymmetric vacimm is obtained by assigning non-trivial 
VEVs to a set of Standard Model singlet fields in the massless string spectrum along F 
and D-flat directions. Some of these fields will appear in the nonrenormalizable terms 
(2.2), leading to effective operators of lower dimension. Their coefficients contain 
factors of order V/M~ 1/10. 

3 Stringent flat directions 

In general, systematic analysis of simultaneously D- and F-flat directions in 
anomalous models is a complicated, non- linear process^. In weakly coupled het- 
erotic string (WCHS) model-building, F-fiatness of a specific VEV direction in the 
low energy efi^ectivc field theory may be proven to a given order by cancellation of 
F-term components, only to be lost a mere one order higher at which cancellation 
is not found. An exception is directions with stringent F-fiatness [5, 6, 20]. Rather 
than allowing cancellation between two or more components in an F-term, stringent 
F-fiatness requires that each possible component in an F-term have zero vacuum 
expectation value. 

When only non-Abelian singlet fields acquire VEVs, stringent fiatness implies that 
two or more singlet fields in a given F-term cannot take on VEVs. For example, in 

^In ref. [26] it is argued that, in addition to flat directions, isolated special points generically 
exist in the VEV parameter space that are not located along flat directions, but for which all D- 
and F-terms are nonetheless zero. Ref. [26] calls upon the proof by Wess and Bagger [27] that 
non- anomalous D-terms do not actually increase the number of constraints for supersymmetric 
flat directions beyond the F-term constraints. However, FI term cancellation, which requires the 
existence of one monomial that is £)-flat for all non-anomalous symmetries, but that carries the 
opposite sign to the Fl-term in the anomalous U{1)a, imposes an additional constraint. Thus, 
without an anomalous U{1), the system of D- and F-equations is not over constraining. In the latter 
case, once a solution to F„i = 0, for all fields s^, is found (to a given order), complexified gauge 
transformations of the fields sji^, that continue to provide a F^ = solution, can be performed that 
simultaneously arrange for non-anomalous D-flatness. Thus, since the F^ = equations impose 
m (non- linear) constraints on rn fields, there should be at least one non- trivial non- anomalous 
-D-flat solution for any set of fields sj^. A parallel proof for non-Abelian field VEVs also exists. 
(Complications to these proofs do arise when different scalar fields possess the same gauge charges.) 
This reduction in apparent total constraints is possible because the F-term equations constrain 
gauge invariant polynomials, which also correspond to non-anomalous F^-flat directions [26] . 



6 



section 4.1, which presents the third through fifth order superpotential for the model 
under consideration, the components of the F-term for <l>45 are (through third order): 



i^*45 = $46*56 + ^-le^se- 



(3.1) 



For stringent F-flatness we require not just that < F^^g >= 0, but that each compo- 
nent within is zero, i.e.. 



Thus, by not allowing cancellation between components in a given F-term, strin- 
gent F-flatness imposes stronger constraints than generic F-flatness, but requires 
significantly less fine-tuning between the VEVs of fields. 

The net effect of all stringent F-constraints on a given superpotential term is 
that at least two fields in the term must not take on VEVs. This condition can be 
relaxed when non-Abelian fields acquire VEVs. Self-cancellation of a single compo- 
nent in a given F-term is possible between various VEVs within a given non-Abelian 
representation. Self-cancellation was discussed in [6] for SU{2) and SO{2n) states. 

A given set of stringent F fiatness constraints are not independent and solutions 
to a set can be expressed in the language of Boolean algebra (logic) and applied 
as constraints to linear combinations of F)-fiat basis directions. The Boolean algebra 
language makes clear that the effect of stringent F-flat constraints is strongest for 
low order superpotential terms and lessens with increasing order. In particular, for 
the model presented herein, stringent fiatness is extremely constraining on VEVs of 
the reduced number of (untwisted) singlet fields appearing in the third through fifth 
order superpotential, in comparison to its constraints on the larger number of singlets 
in the model of [7]. 

One might imagine that stringent F-flatness constraints requires order-by-order 
testing of superpotential terms. This is, in fact, not necessary. All-order stringent 
F-fiatness can actually be proven or disproven by examining only a small finite set 
of possible dangerous (i.e., F-fiatness breaking) superpotential terms. Through a 
process such as matrix singular value decomposition (SVD)^, a finite set of superpo- 
tential terms can be constructed that generates all possible dangerous superpotential 
terms for a specific D-fiat direction. This basis of gauge-invariants can always be 
formed with particular attributes: (1) each basis element term contains at most one 
unVEVed field (since to threaten F-fiatness, a gauge-invariant term, necessarily with- 
out anomalous charge, can contain no more than one unVEVed field); (2) there is at 
most one basis term for each unVEVed field in the model; and (3) when an unVEVed 
field appears in a basis term, it appears only to the first power. The SVD process 
generated a possibly threading basis of superpotential terms for several models (see 
for example [5, 6, 9, 17, 29, 30, 31]). 

SVD fortran subroutine is provided in [28]. 



< $46*56 >= 0> < *46*56 >= 0. 



(3.2) 
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To appear in a string-based superpotential, a gauge invariant term must also follow 
Ramond-Neveu-Scliwarz worldslieet charge conservation rules. For free fermionic 
models these rules were generalized from finite order in [24, 32] to all-order in [17]. 
The generic all order rules can be applied to systematically determine if any product 
of SVD-generated F-fiatness threatening superpotential basis elements survive in the 
corresponding string-generated superpotential. If none survive, then F-flatness is 
proven to all finite order. This technique has been used to prove F-flatness to all 
finite order for various directions in several models [5, 6, 9, 17, 29, 30, 31]. Alternately, 
if any terms do survive, the lowest order is determined at which stringent F-fiatness 
is broken. 

How should stringent (especially all-order) flat directions be interpreted in com- 
parison to general (perhaps flnite order) flat directions? All-order stringent flat direc- 
tions contain a minimum number of VEVs and appear in models as the roots of more 
fine-tuned (generally finite-order) flat directions that require speciflc cancellations 
between F-term components. The latter may involve cancellations between sets of 
components of different orders in the superpotential. 

All-order stringent flat directions have indeed been discovered to be such roots 
in all prior free fermionic heterotic models for which we have performed systematic 
flat direction classiflcations. However, the model presented herein appears to lack 
any stringent flat directions, at least within the expected range of VEV parameter 
space. We have reached this conclusion after employing our standard systematic 
methodology for D- and F-flat direction analysis. 

4 The String Model 

In this section we give the details of our model. The boundary condition ba- 
sis vectors beyond the NAHE-set and the one-loop GSO projection coefficients are 
shown in eq (4.1) and eq. (4.2), respectively. 
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With the choice of generalised GSO coefficients: 
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(4.2) 



In matrix (4.2) only the entries above the diagonal are independent, while those 
below and on the diagonal are fixed by the modular invariance constraints. Blank 
lines are inserted to emphasise the division of the free phases between the different 
sectors of the reahstic free fermionic models. Thus, the first two lines involve only the 
GSO phases of c(^^^f^). The set {1, S} generates the = 4 model with S being the 

space-time super symmetry generator and therefore the phases c(^^^ are those that 
control the space-time supersymmetry in the supcrstring models. Similarly, in the 
free fermionic models, sectors with periodic and anti-periodic boundary conditions, 
of the form of 6j, produce the chiral generations. The phases c(^^'J determine the 
chirality of the states from these sectors. 

Both the basis vectors a and (3 break the S'O(IO) symmetry to S0{&) x SO {A) 
and the basis vector 7 breaks it further to 5f/(3) x U{l)c x SU{2) x U{1)l. The 
basis vector a is symmetric with respect to the sector hi and asymmetric with respect 
to the sectors 62 and 63, whereas the basis vector j3 is symmetric with respect to 62 
and asymmetric with respect to hi and 63. As a consequence of these assignments 
and of the string doublet-triplet sphtting mechanism [33], both the untwisted Higgs 
colour triplets and electroweak doublets, with leading coupling to the matter states 
from the sectors hi and 62, are projected out by the generalised GSO projections. 
At the same time the untwisted colour Higgs triplets that couple at leading order to 
the states from the sector 63 are projected out, whereas the untwisted electroweak 
Higgs doublets remain in the massless spectrum. Due to the asymmetric boundary 
conditions in the sector 7 with respect to the sector 63, the leading Yukawa couphng 
is that of the up-type quark from the sector 63 to the untwisted electroweak Higgs 
doublet [4]. Hence, the leading Yukawa term is that of the top quark and only its 
mass is characterised by the electroweak VEV [4]. The lighter quarks and leptons 
couple to the light Higgs doublet through higher order nonrenormalizable operators 
that become effective renormalizable operators by the VEVs that are used to cancel 
the anomalous U{1)a -D-term equation [4]. The novelty in the construction of ref. [7], 
and in the model of eq. (4.1), is that the reduction of the untwisted Higgs spectrum is 
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obtained by the choice of the boundary condition basis vectors in eq. (4.1), whereas 
in previous models it was obtained by the choice of fiat directions and analysis of the 
superpotential [10]. 

The final gauge group of the string model arises as follows: in the observable 
sector the NS boundary conditions produce gauge group generators for 

SU{3)c X SU{2)l X U{l)c x ^7(1)^ x t/(l)i,2,3 x t/(l)4,5,6 . (4.3) 
Thus, the ,50(10) symmetry is broken to SU{3) x SU{2)l x U{1)c x U{1)l, where, 

C/(l)c = IVC/(3)c Qc^Y^Qm. (4.4) 

j=i 

5 

U{1)l = IVC/(2)i QL^Y^Qi^')- (4-5) 

The flavour 5*0(6)^ symmetries are broken to [/(l)^"*"" with (n = 0, ■•■,6). The 
first three, denoted by U{l)j {j = 1,2,3), arise from the world-sheet currents fj^fj^* . 
These three f/(l) symmetries are present in all the three generation free fermionic 
models which use the NAHE set. Additional horizontal U{1) symmetries, denoted 
by U{l)j {j — 4,5,...), arise by pairing two real fermions from the sets {y^''"'^}, 
{y^'^ , u)^'^} and {lD^''"'^}. The final observable gauge group depends on the number 
of such pairings. In this model there are the pairings y^y^, y^uT' and (D^tt)^, which 
generate three additional [/(I) symmetries, denoted by f/(l)4,5,6- 

It is important to note that the existence of these three additional U{y) currents 
is correlated with the assignment of asymmetric boundary conditions with respect 
to the set of internal world-sheet fermions {y, a;||/, d)}^''"'^, in the basis vectors that 
extend the NAHE-set, {a. /?, 7}. This assignment of asymmetric boundary conditions 
in the basis vector that breaks the 50(10) symmetry to 5*0(6) x 5*0(4) results in 
the projection of the untwisted Higgs colour-triplet fields and preservation of the 
corresponding electroweak-doublet Higgs representations [33]. 

In the hidden sector, which arises from the complex world-sheet fermions ^^'"^ ■, 
the NS boundary conditions produce the generators of 

^f/(2)i,2,3,4 X SU{^)h, X U{\)h, ■ (4.6) 
U{\)ii^ corresponds to the combinations of the world-sheet charges 

Qh.^YQW). (4.7) 

The model contains several additional sectors that may a priori produce space- 
time vector bosons and enhance the gauge symmetry, which include the sectors C = 
1 + 61 + 62 + &3 and 1 + 5' + q; + /3 + 7. Additional space-time vector bosons from 
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these sectors would enhance the gauge symmetry that arise from the space-time 
vector bosons produced in the Neveu-Schwarz sector. However, with the choice of 
generahsed GSO projection coefficients given in eq. (4.2) all of the extra gauge bosons 
from these sectors are projected out and the four dimensional gauge group is given 
by eqs. (4.3) and (4.6). 

In addition to the graviton, dilaton, antisymmetric sector and spin-1 gauge 
bosons, the Neveu-Schwarz sector gives one pair of electroweak Higgs doublets 
and h^; six pairs of 5*0(10) singlets, which are charged with respect to f/(l)4,5,6; 
three singlets of the entire four dimensional gauge group. A notable difference as 
compared to models with unreduced untwisted Higgs spectrum, like the model of ref. 
[3], is that the 5*0(10) singlet fields, which are charged under f/(l)i,2,3, are projected 
out from the massless spectrum. The three generations are obtained from the sec- 
tors bi, 1)2 and 63. The model contain states that are vector-like with respect to the 
Standard Model and all non-Abelian group factors, but may be chiral with respect 
to the [/(I) symmetries that are orthogonal to the 5*0(10) group. The full massless 
spectrum of the model is detailed in Table 1 at the end of this paper. 

As a final note we remark that the boundary conditions with respect to the 
internal world-sheet fermions of the set {y,uj\y,uo}^'"''^ in the basis vectors a, j3 and 
7, that extend the NAHE-set, are similar to those in the basis vectors that generate 
the string model of ref. [3] , with the replacements 



The world-sheet fermions {y, cul^, cD}^' "'^ correspond to the compactified dimen- 
sions in a corresponding bosonic formulation. The substitutions in (4.8) are aug- 
mented with suitable modifications of the boundary conditions of the world-sheet 
fermions j?/;^' "'^, f;^' "'^, 0^'' "'®}, which correspond to the gauge degrees of freedom. 
The effect of these additional modifications is to alter the hidden sector gauge group. 
While the substitutions in (4.8) look innocuous enough, they in fact produce sub- 
stantial changes in the massless spectrum and, as a consequence, in the physical 
characteristics of the models. With regard to the flat directions of the superpoten- 
tial, the effect of these changes on the untwisted states will be particularly noted. 

4.1 Third through Fifth Order Superpotential 

The three singlets of the entire four dimensional gauge group are obtained from: 



«(y¥) 



" 7(yV). 



(4.8) 



6 
6 



X^^^cu^cu^lO > 



We show below the cubic through quintic order superpotential terms. 
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Trilinear superpotential: 

W3 = N'^L,,h + ulQ^h + H^Hrh + H^Hrh + 

+ ^2{H2H2 + HiqHiq + HiiHii) 

+ ^siHsHs + H4H4 + + HqHq + HiHi) 

+ ^3($f l»f + $f l>f ) 

+ $45($46l>56 + ^le^Se) + ^45(^46^56 + ^46^56) 
+ $l5(l'46^56 + ^'aA) + ^45(^46^56 + ^'aA) 

+ $'45 {{^ff + {^ff) + )^ + im^f) 

+ ^45-^12-^^13 + $46-f^l4-f^l5 + ^'^qHiqHij 

+ ^fHgHn + ^^{HA + HsHio) + H1H4H10 + H^H^Hs (4.9) 
Quartic superpotential: 

Wa = QiUiH4H5 + Q2U2H4He + LiN^H4H5 + L2N^H4HQ (4.10) 
Quintic superpotential: 

W5 = QiHsL,H,^2 + Q2H^L2H^^i + Q^ulH^H^H^^ + Qsu^H^HrHs 
+ dlu1H,H,^2 + d1H,H,^,,V2 + d'^^ulH^H^i, + d',HsH,^',,V, 
+ H-iH4^HiH^HiQ + HSH4H2H3HS + H^HiH2H2,^2'^ + H^H^^"^ HnHg 
+ HsHs^2^ HgHio + LsHiN^HjHio + L^H2N^H'jH^ + H4^H4^^qHsHs 
+ H,H,<^4,N'[V2 + i/4i^45';6^2K5 + H4H4^',,H,oH,o + H^H^H^H^H^, 
+ H4H4H4H2HS + -ff4-/?4-/?ii5^2*^2^ + H^H^^i^ HiiHg -f- H/^H^^'^^ H^HiQ 
+ H4HiC^2H03Hs + H4H2^^Q^2^ Hio + HaHa^^qHsHs + HaHa^^qHiqHio 
+ H4H1H1H1H1Q + H4^HiHiH2Hg + H4^HiH2H2HiQ + H^HiHjHjHif) 
+ H4^HiHqHqHio + Hi^HiH^H^HiQ + H4HiC,2'^2^ Hs + HaHi^i'^'^i'^Hio 
+ H4Hi^2^^2'^Hio + HaHiHuHioHh + H4H1H10H10H10 + H^HiHioHgHg 

+ HaHiHiqHsHs + H^HiHiqHiqHiy + H4H2H2H2H8 + -^4-^2*J'56*^2'^-^10 

+ H4H2H7HiHs + H4H2HQHQHS + -^4-f^2-f^5-f^5-f^8 + -f^4-f^2*^'l ^-^^8 

+ H4H2^2^ ^2^ Hs + H4H2H11HSH11 + H4H2H1QHSH10 + H4H2H9HSH9 

+ H4H2H8H8HS + HiHiH^H2^f + i?ii?ii?2i?2<^>l5 + i?iiyi$;6<i>f $f 

+ ^i^i^le^f + HiHi<^fHuHg + H^H^^fHsHw + HiH2H2H2^f 

+ H,H2<^'^,HsHw + HiH2H^H^^f + HiH2HeHe^f + H,H2H,H,^f 

+ HiH2^f^f^f + HiH2^f^f^f + HiH2^f^f^f + HiH2^fHnHi 

+ HiH2^fHioHio + HiH2^fH9Hg + HiH2^fHsHs + H^H^H2H2^'4^ 
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+ ^2^2$56^f + H2H2^fHnHg + H2H2^fHsHio + H2H2^'^f,^f 

+ H2H2^'^^^f^f + H2H2^',eHi2His + ^'^^^f^fHi^Hn + ^'^^^f^fH.eHir 

+ <^,,N',VsH^,Hn + $ 45iVrV^2i^io^io + ^seiVs^gi^ii^ii + ^,eN^VsH,oH,o 

+ ^,,^f^fH,,H,, + ^,,^f^fH,,H,, + N^V,^fH,oHs + N^^fHuHsVu 

+ HjHj^"^ HuHq + HtHj^2'^ H^HiQ + HqHq^i'^ HuHq + HqHq^2^ HsHio 

+ HuHi2HqHi3 + ^2^ HiiH^HiiHiQ + ^2^ HiqHsHiqHw + ^2^ HgHsHioHg 

+ ^fSsHsH.oHs (4.11) 



5 Flat directions 

The model possesses nine local U{1) symmetries, eight in the observable part and 
one in the hidden part. Six of these are anomalous: 

TrC/i = TrC/2 = -TrC/3 = 2TrC/4 = -2TrC/5 = 2TrC/6 = -24. (5.1) 

U{1)l and ^/(l)c of the 50(10) subgroup are anomaly free. Consequently, the weak 
hypercharge and the orthogonal combination, U{l)z', are anomaly free. The hidden 
sector U{1)hi is also anomaly free. 

Of the six anomalous U{l)s, five can be rotated by an orthogonal transforma- 
tion to become anomaly free. The unique combination that remains anomalous is: 
Ua = kJ2j[T^^U{l)j]U{l)j, where j runs over all the anomalous U{l)s and /c is a nor- 
malisation constant. For convenience, we take k — ^ and therefore the anomalous 
combination is given by: 

Ua = -2Ui - 2U2 + 2f/3 - f/4 + f/5 - Ua, TtQa = 180. (5.2) 

We note that the anomalous U{1)a combination in (5.2) is similar to that of the 
model of ref. [3], aside from the changes of the sign of the traces over U{l)i2 and 

^(1)5- 

The five rotated non-anomalous orthogonal combinations are not unique, with 
different choices related by orthogonal transformations. One choice is given by: 

U\ ^ U1-U2 , U'2 = U1 + U2 + 2U3, (5.3) 
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C/'3 ^ U4 + U5 , U\ = C/4 - C/5 - 2C/6, 
C/'s = C/i + C/2 - C/3 - 2C/4 + 2C/5 - 2C/6. 



(5.4) 
(5.5) 



Thus, after this rotation there are a total of eight U{l)s free from gauge and gravi- 
tational anomalies. 

A basis set of (norm-squares of) VEVs of scalar fields satisfying the non-anomalous 
D-fiatness constraints (1.4) can be created en masse [15, 16, 17]. The basis directions 
can have positive, negative, or zero anomalous charge. In the maximally orthogo- 
nal basis used in the singular value decomposition approach of [16, 17], each basis 
direction is uniquely identified with a particular VEV. That is, although each basis 
direction generally contains many VEVs, each basis direction contains at least one 
particular VEV that only appears in it. 

A physical D-flat direction i^phys, with anomalous charge of sign opposite that of 
the FI term is formed from linear combinations of the basis directions, 

# basis dirs. 

-Dphys = XI ^i^i^ (^-6) 
i=l 

where the integer coefficients are normalised to have no non-trivial common factor. 

In our notation, a physical flat direction (5.6) may have a negative norm-square 
for a vector-like field. This denotes that it is the oppositely charged vector-partner 
field that acquires the VEV, rather than the field. Basis directions themselves may 
have vector-like partner directions if all associated fields are vector-hke. On the 
other hand, if in particular, the field generating the VEV uniquely associated with a 
basis direction docs not have a vector-like partner, that basis direction cannot have 
a vector-like partner direction. 

In pursuit of physical all-order fiat directions for this model, we first examined 
directions formed solely from the VEVs of non-Abehan singlet fields. An associated 
maximally orthogonal basis set, denoted by {2^j=i to 13}) containing only non-Abelian 
singlet VEVs is shown in Table 2. a. The respective unique VEV fields of these basis 
directions are identified in Table 2.b. Examination of Tables 2. a and 2.b reveals that 
no physical D-fiat directions can be formed solely from VEVs of non-Abelian singlet 
fields. Since the FI term ^ (1.4) is positive for this model, with TrQ^ = 180, a 
physical flat direction must carry a negative anomalous charge. However, of the 13 
singlet Z^-flat basis directions, three carry anomalous charge of +15, +30, +30 while 
the remaining ten do not carry anomalous charge. Further, the unique VEVed fields 
for the 3 basis directions with positive anomalous charge do not have corresponding 
vector-like partner fields. Hence, there are no vector-like paired basis directions with 
negative anomalous charge. Thus, Tables 2. a and 2.b imply that one or more fields 
carrying non-Abelian charges must also acquire VEVs in physical D-flat directions. 
This result is, in itself, not necessarily unexpected, as non-Abelian VEVs have been 
required for physical (all-order) fiat directions in other quasi-realistic free fermionic 
heterotic models in the past, for example [9]. 



14 



Thus, wc expanded our flat direction search to include VEVs of both non-Abelian 
singlet fields and non-Abelian charged fields. Our chosen set of 50 maximally or- 
thogonal D-flat basis directions for both non-Abelian singlet VEVs and non-Abelian 
charged VEVs, denoted by {I^i=i to 50}, is presented in Table 3. a. The respective 
unique field VEVs identified with these basis directions are given in Table 3.b. In 
this enlarged basis the anomalous charges are given in units of (^5-) and the direc- 
tions containing only singlet VEVs are rotations of those in Table 2. a. 

Nine of the 50 directions, denoted -Di=i,...,9, carry one or two units of nega- 
tive anomalous charge. Twenty basis directions, denoted Diq through D29, carry 
no anomalous charge. Twenty-one basis directions, denoted D30 through D^q, 
carry one or two units of positive anomalous charge. All basis directions possess- 
ing negative anomalous charge contain SU{3)c ® SU{2)l charges or hidden sector 
SU (4) ® 111=1 charges. (Thus, this basis set also reveals that anomaly cancel- 

lation will necessarily break one or more non-Abelian local symmetries.) All of the 
$ fields, the Hi to 11 fields and h have vector-hke pairs. Thus, physical flat directions 
can have negative components for any of these. A subset of these fields, specifically 
^46, ^45, ^56) -f^4,5,6,7, has VEVs appearing in multiple basis directions. The only 
non- vector-like field with a VEV that appears in multiple directions is Cg. 

DiQ through D^-j and D22 are composed solely of varying combinations of the 
vector-like fields. Hence, all of these basis directions have corresponding vector-like 
partner basis directions, Di = —D^, for which the VEV of each field is replaced by the 
VEV of the vector- like partner field. Thus, in a physical fiat direction (5.6), each of 
the respective integer coefficients aio through ayj and 022, may be negative, positive, 
or zero. 

Note that D-j, Dg, Dg and D20 are vector-like except for their Cg components. 
Thus, each of ay, ag, Qq and 020 may be negative, positive, or zero in a physical D-Hat 
direction, so long as the net norm-square VEV of Cg is non-negative.^ The remaining 
basis directions contain at least one unique non-vector-like field VEV. Thus, in a 
physical fiat direction, the coefficients of the remaining basis directions must be non- 
negative. 

What does this mean for a physical D-Hat direction formed as a linear combination 
of the basis directions? For a physical flat direction there are, thus, two speciflc 

constraints on the coefficients and one general set of non-negative norm-square 
constraints on a subset of the a.;. First, negative anomalous charge for a fiat direction 

^Note that non- vector-like fields, such as e§, that appear in multiple directions with some basis 
directions having positive and some having negative norm-square components, are common in this 
process. Further, some models explored in the past have had (at least) one basis direction with 
two (or more) field VEVs unique to it and with norm-square VEVs with difi'ering signs. This latter 
type of basis direction can never appear in a physical direction and, hence, implies that the fields 
unique to it can never appear in a _D-flat direction. (If all of the norm-squares of the fields unique 
to a basis direction were initially negative, then these signs, along with those of the norm-squares of 
any vector-like field VEVs in that basis direction, could all be changed together to allow the basis 
direction to appear in a physical direction.) 
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requires 

2 9 44 50 

-2E«'* + E «i + 2 E «i < 0- (5-7) 

t=l j=3 i=30 j=45 

Second, a non-negative norm-square VEV for requires 

2 6 9 19 

-6 E - 3a3 - 6 E - 2a7 - 6 E - 2 E a, - + «2i 

1=1 i=4 i=8 i=18 

24 

-2 E Oj - a25 - 2a26 + 2a27 - 2a28 + 2a29 + 6030 + 6032 + a^s 

i=23 

40 47 49 

+3^ ai + 6a42 -^6^0^-^2^01-1- 6059 > 0. (5.8) 

1=39 i=45 i=48 

Last, for the set of non- vector- like fields that are each identified with a respective 
unique D-Hat direction, the general set of non-negative norm-square VEV constraints 
is 

Ui > for i = 1 to 6, 18, 19, 21, 23 to 50. (5.9) 

At low orders, each individual supcrpotential term also induces several stringent 
F-term constraints on the a, coefficients of physical fiat directions. As stated prior, 
the set of constraints from superpotential terms with only singlet fields translate 
into the requirement that two or more singlet fields in a given superpotential term 
cannot take on VEVs. For the model under investigation, constraints from third order 
superpotential terms are especially severe. For this model, all six $ singlet fields and 
their vector-like partners appear in third order superpotential terms (specifically, the 
sixth and seventh fines) of (4.9). Stringent F-fiatness from these terms forbids at 
least 8 of the 12 singlet fields from acquiring VEVs. 

For example, when solely third order stringent F-flatncss constraints are applied 
to the six pairs of $ vector-like singlets (and no F-flatness constraints are applied to 
the non-Abelian states), there are just nine solution classes that allow the maximum 
of 4 singlet VEVs. (Flat directions in any of these nine classes are defined by their 
respective non-Abefian VEVs.) 

For three of these nine singlet third order flatness classes, the VEVs are of two 
fields and their respective vector-like partners: either, 

< $45 >, < ^^5 >, < ^45 >, < ^45 > 7^ 0, or (5.10) 

< l'46 >, < $46 >, < $46 >, < ^46 > 7^ 0> Or (5.11) 

< l>'56 >, < $56 >, < ^56 >> < ^56 > 7^ 0. (5.12) 

Higher order stringent flatness constraints can further reduce the allowed number 
of singlet VEVs of each of these solutions. Further, a component of a D-Hai basis 
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direction in Table 3. a only specifies the difference between the norm-squares of the 
VEV of a given field and of the given vector-like partner field (if it exists). Completely 
chargeless VEVs solely involving a field $i and its vector-like partner $j such that 
|<$i>P = |<l>i>P can always be added to a physical D-flat direction. 
However, it is preferable for higher order F-fiatness to impose that a field and its 
vector-partner do not simultaneously acquire VEVs. Hence, these three solutions 
effectively allow only two unique singlet fields to acquire VEVs. 

The next three classes of singlet solutions do allow up to four distinct singlet fields 
to acquire VEVs: either. 



For the three remaining solution classes, the fields in (5.13), (5.14) and (5.15), 
are respectively replaced with their vector-like partner fields. For any of these nine 
stringent F-flat choices, no other $ singlet fields can acquire VEVs. 

Any of the constraints on allowed and disallowed VEVs, such as the above, can 
be re-expressed in terms of constraints on the Oj coefficients specifying the basis 
directions contributions to a physical D-fiat direction. For example, setting < $45 >= 
would require 



To systematically investigate physical D-fiat directions with non-Abelian VEVs, 
over a course of eight months we generated and examined physical D-flat directions 
composed of from 1 to 6 basis directions. Under the assumption that all VEVs 
of physical flat directions are nearly of the same order of magnitude, we allowed 
coefficients of to 20 for the non-vector-like basis directions and coefficients of -20 
to 20 for the vector-like basis directions. 

To be classified as a physical D-Hat direction, a linear combinations of basis direc- 
tions needed to obey (5.7-5.9) and was, of course, also required to have non-Abelian 
D-flatness. (The general process by which we enforced non-Abelian D-flatness fol- 
lowed that presented in [16, 17].) Each resulting physical D-flat direction was then 
tested for stringent F-flatness from all third order through fifth order superpotential 
terms and additionally for some key sixth order superpotential terms. ^ 

^While only the third through fifth order superpotential is given in section 4.1, we have generated 
the complete superpotential to eighth order and can generate it to any required order. 



< $45 >, < >, < $46 >, < $l6 >^ 0, or, 

< $45 >, < $'56 >, < >< $^5 >^ 0, or, 

< $46 >, < $56 >, < $56 >> < "^16 >^ 0- 



(5.13) 
(5.14) 
(5.15) 



4 

4:ai + a2 + 2^ai + 805 + 2aQ + — as — ag + Oio + Oie 

•i=3 

— ai8 -|- 2ai9 + 0.20 ~ 0,21 + 0-23 ~ 2a27 ~ Ct28 + 0.29 + 4a3o 

35 37 43 

+a3i - 032 ^ - 2 ^ - 2039 + a4o-2^ai + 044 

j=33 i=36 i=41 

— 4ci45 + 2a,4g — 0,47 — 3049 — 050 = . 



(5.16) 
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Following the SVD method discussed earlier in section 3 and described in [5, 6, 20], 
we had planned to then test for possible all-order stringent F-fiatness, the subset of 
physical D-Aat directions that had proven stringently F-fiat to at least fifth or sixth 
order. Based on all of the prior models we had investigated, we had expected to 
find around four to six physical D-Hsit directions that were, in fact, stringently F- 
flat to all finite order. However, in contrast we discovered that no physical D-fiat 
directions that we had generated even kept stringent F-fiatness through sixth order. 
So there were no physical D-Hat directions to examine for all-order testing! For this 
model, with its reduced set of singlet fields from the untwisted sector, not even self- 
cancellation of non-Abehan terms could provide stringent F-flatness through sixth 
order for any of these physical D-flat directions. 

We will continue for several years a search for F-fiatness past sixth order for 
physical D-flat directions in this model that are comprised of seven or more basis 
directions. After January 2008, the search will be conducted on a 128-node com- 
puter system of quad-processors to quicken the pace of the investigation. However, 
a continued null result is likely: since each of our basis directions contains a unique 
field VEV, increasing the number of non-zero aj coefficients linearly increases the 
minimum number of unique field VEVs. With each increase in number of basis di- 
rections composing a physical F)-fiat direction, the probability of obtaining stringent 
F-flatness much beyond sixth order further decreases. 

6 Conclusions 

The string models in the free fermionic formulation gave rise to a large class of 
quasi-realistic string models, including three generation models that produce solely 
the MSSM spectrum in the observable Standard Model charged sector of the effective 
low energy fleld theory. As such the free fermionic models provides the arena to 
study how string theory may be related to the observed particle data. In turn, the 
properties of the models that make them attractive from the point of view of the 
phenomenological data may be instrumental to uncover unexpected properties of 
string theory. 

In this paper we stumbled upon such a possible novel string feature with regard 
to supersymmetry breaking. A general expectation from super symmetric quantum 
field theories is that if supersymmetry is unbroken at the lowest order in perturbation 
theory, then it cannot be broken at higher perturbativc orders in quantum perturba- 
tion theory. The model that we presented in this paper opens up the possibility that 
string theory may afford other options. 

The model is a quasi-realistic three generation string model in the free fermionic 
formulation that shares many of the characteristics of previous quasi-realistic free 
fermionic models. It contains three chiral generations charged under the standard- 
like model subgroup of the underlying 5*0(10) symmetry of the NAHE-set. The 
weak-hypercharge possesses the canonical GUT embedding and the model predicts 
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sin^ = 3/8 at the unification scale. The Higgs spectrum contains one untwisted 
electroweak doublet pair that couples at the cubic level of the superpotential to the 
top quark. Like numerous other quasi-realistic string models, the model contains an 
anomalous U{1) symmetry which generates a Fayet-Iliopoulos L>-term. The Fayet- 
Iliopoulos term would, in general, break supersymmetry, unless there is a direction 
in the scalar potential which is both D- and F-flat. If such a direction exists it 
will acquire a VEV, cancel the Fayet-Iliopoulos D-term and restore supersymmetry. 
The general expectation, due to the corresponding results in supcrsymmctric point 
quantum field theories, is that a supersymmetric vacuum does exist due to the fact 
that string spectrum is Bose-Fermi degenerate and possesses N — 1 supersymmetry. 

Indeed, in all the previously studied quasi-realistic free fermionic models such 
supersymmetric flat directions were found. Moreover, all previous models yielded 
the so called stringent flat directions that can be shown to be exact, i.e. flat to all 
orders of nonrenormalizable terms. The distinct feature of the model discussed in this 
paper is that it does not admit such stringent flat directions. In this model no physical 
D-Hat direction that we generated kept F-flatness through sixth order. We speculate 
that only stringent flat directions can be flat to all orders of nonrenormalizable terms. 
If it is validated, this would indicate that this model, therefore, appears to have no 
D-flat directions that can be proven to be F-flat to all orders, other than by order 
by order analysis. 

As we discussed this outcome may be a general result of the assignment of bound- 
ary conditions to the internal world sheet fermions, which results in the projection 
of two of the untwisted Higgs pairs. If a non-vanishing F-tcrm docs exist, the 
implication would be that supersymmetry remains unbroken at flnite order. The 
Fayet-Ihopoulos term that breaks supersymmetry is generated at the one-loop level 
in the perturbative string expansion. On the other hand the string spectrum is Bose- 
Fermi degenerate and possesses = 1 space-time supersymmetry at the classical 
level. This would suggest that, contrary to the expectation from supersymmetric 
point quantum fleld theories, perturbative supersymmetry breaking in string theory 
may ensue. Furthermore, the modular invariant one-loop partition function vanishes 
and, hence, the cosmological constant vanishes at the one loop level as well. The 
model presented here may therefore represent an example of a quasi-realistic string 
vacuum with vanishing one-loop cosmological constant and perturbatively broken su- 
persymmetry. Furthermore, the asymmetric boundary conditions of the string model 
given in eq. (4.1) project all the geometrical moduli of the underlying Z2 x Z2 orbifold 
[34, 35]. Absence of supersymmetric flat solutions would imply that the supersym- 
metric moduli are flxed as well in this model. Examining the hidden sector gauge 
group of the model, given in eq. (4.6), we note that it contains SU{2)^ and satisfles 
the conditions for the dilaton race-track stabilisation mechanism [36]. It remains 
without saying that many issues still need to be further explored and understood to 
ascertain the claims of this paper. Nevertheless, the free fermionic models continue 
to generate intriguing and exciting results that may at the end prove to be relevant 
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to the elucidation of the connection between string theory and the real world. 
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Table 2. a. D-Flat direction basis of non-Abelian singlet fields. 



Column 2 specifies the anomalous charge and columns 3 through 16 specify the 
norm-square VEV components of each basis direction. The six fields and i^5,6,7 
carry hypercharge, the remaining do not. A negative component indicates the vector 
partner of a field (if it exists) must take on VEV rather than the field. 
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Table 2.b. Unique VEV associated with each non-Abehan singlet field D-Flat basis 
direction. 
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Table 3. a D-Flat direction basis of all fields. 
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Table 3. a continued. 
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Table 3. a continued. 
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Table 3.b. Unique VEV associated with each D-Flat basis direction. 
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